Quantum correlations in the soliton collisions 
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We study quantum correlations and quantum noise in the soliton collision described by a general 
two-soliton solution of the nonlinear Schrodinger equation, by using the back-propagation method. 
Our results include the standard case of a sec/i-shaped initial pulse analyzed earlier. We reveal 
that double-hump initial pulses can get more squeezed, and the squeezing ratio enhancement is due 
to the long collision period in which the pulses are more stationary. These results offer promising 
possibilities of using higher-order solitons to generate strongly squeezed states for the quantum 
information process and quantum computation. 
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I. INTRODUCTION 

Light squeezing is an important physical concept that 
continues to attract attention of researchers due to its 
potential for implementing quantum information pro- 
cessing and quantum computing. As an alternative to 
single-photon schemes, the demonstration of the Ein- 
stein, Podolsky, and Rosen (EPR) paradox and quan- 
tum teleportation with continuous variables have been 
realized experimentally by using the entanglement from 
squeezed states 0,0- Moreover, experimental progress 
in the study of various quantum information processing 
with squeezed states ggnerated from optical fibers has 
recently been reported |^ 0,0,1^. To increase the entan- 
glement fidelity of continuous variables, the enhancement 
of squeezing effect becomes very vital. 

Original proposals to generate squeezed states from op- 
tical fibers are based on the use of the fundamental soli- 
tons supported by the Kerr nonlinearity of silica glass. 
Temporal (pulse) solitons in optical fibers are described 
by the nonlinear Schrodinger equation (NLSE) that can 
exhibit quadrature-field squeezing 0, Hi |9|, [l^, El as well 
as amplitude squeezing FalT^. and both intra-pulse and 
inter-pulse correlations |l4| . Besides the exact quantum 
soliton NLSE solutions constructed by using the Bethe 
ansatz p"(l | , the quantum properties of temporal solitons 
are well described by the linearization approach 11] for 
the average photon number as high as 10^. Based on this 
linearization approach, many different numerical meth- 
ods have been developed during the past two decades in 
order to study quantum noise associated with nonlinear 
pulse propagation, including the positive-P representa- 
tion 7, 8], back-propagation method [lHI |. and the cumu- 
lant expansion technique |l6j |. 

Experimentally, the soliton squeezing from a Sagnac 
fiber interferometer, 1.7 dB below the shot noise, was first 
observed in 1991 by Rosenbluh and Shelby Since 
that, larger quadrature squeezing from fibers has been 
obtained with a gigahertz Erbium-doped fiber lasers that 
allow to suppress the guided acoustic-wave Brillouin scat- 



tering, and 6.1 dB noise reduction below the shot noise 
has been reported attempt to enhance the 

soliton squeezing effect, one may increase the energy of 
an optical soliton enhancing the importance of nonlinear 
effects and, employing this idea, 7.1 dB photon-number 
squeezing has been demonstrated by using spectral fil- 
ters [Tgt . 

However, it is known that the basic model of the pulse 
propagation in optical fibers described by NLSE pos- 
sesses more general iV— soliton solutions which can be 
obtained, for example, by applying the inverse scattering 
transform 20] . As was demonstrated, such higher-order 
{N = 2,3,...) solitons can be more squeezed since they 
contain times of the energy than the fundamental 
soliton [T^l2lll2^ : as an example, up to 8.4 dB enhance- 
ment was predicted for the N = 2 soliton states j23|. 

It is important to mention that all previous studies 
of the soliton quantum noise and quantum squeezing of 
higher-order solitons have employed a very special case 
of two-soliton states generated by a simple sech-like in- 
put pulse. However, a general NLSE solution describing 
the 7V-soliton state is characterized by N free parame- 
ters which can be controlled independently. In this pa- 
per, we develop the theory of quantum noise and quan- 
tum squeezing in the context of the multi-soliton states 
and apply it to study squeezing of the general A^-soliton 
bound states of NLSE. In particular, we reveal that 
the conventional sech-like single-hump pulses are not the 
most suitable pulses for generating highly squeezed states 
and, using the case of the general two-parameter solution 
for N — 2 solitons as an example, we show that an in- 
put double-hump soliton is better for generating strongly 
squeezed states, even such solitons have the same en- 
ergy as the single-hump pulses. The enhancement of the 
squeezing effect is explained by the long collision period 
of a double-hump soliton, consequently the pulse profile is 
more stationary for getting squeezed. Since these double- 
hump solitons have also been generated in the fiber laser 
systems, more strongly squeezed state from optical fibers 
are expected to be realized with the current technology. 
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FIG. 1: The contour plots the evolution of the N = 2 solitons 
with (a) r;i : 772 = 1 : 3, and (b) 771 : 7^2 = 1 : 2. The 
insets show the initial soliton profiles at z = 0. In all cases, 
t;i + 77 = 2. The straight lines mark the propagation distances 
shown in Fig. |21 



II. TWO-SOLITON BOUND STATES 



When 771 : 772 = 1 : 3, the soliton solution at z 
specific, sech-like single-hump initial profile, 

C/(0,i) 



has a 



sech(i) ' 

as shown in the insert of Fig. ^a). However, when the 
ratio of 771/772 becomes larger than 1/3, the initial profile 
of the N = 2 soliton solution becomes double-humped, 
as shown in the insert of Fig.^b) for the special case of 
?yi : ??2 = 1 : 2. In spite of such a difference in the soliton 
profiles, the soliton energy defined as 

P = J \U{z,t)\dt^, 

remains the same for the full set of the N — 2 soliton 
solutions, i.e. P = 8 for 771 -|- 772 = 2 and arbitrary ratio 



To describe the pulses propagating in optical fibers 
with the anomalous dispersion and Kerr-type nonlinear- 
ity, one employ the NLSE model written for the normal- 
ized variables, z and t, 



.dU{z,t) d^U{z,t) 



dz 



9t2 



+ \U{z,t)\'^U{z,t)^0, (1) 



where U{z, t) is the pulse envelope. According to the re- 
sults of the inverse scattering transform [23, this equa- 
tion possesses exact solutions describing interaction of N 
solitons, which are characterized by a set of complex vari- 
ables {Aj, Cj}, j = 1, 2, . . . , TV. The complex Aj and Cj 
are the "poles" and the "residues" of the corresponding 
scattering data [23|. In particular, such solutions de- 
scribe the so-called TV— soliton bound states (also called 
"breathers") when all solitons have vanishing velocities 
at the infinity, and their interaction leads to the forma- 
tion of a spatially localized but time-periodic state. The 
full set of such solutions for the N = 2 soliton bound 
states can be written as |2J, ^91 , 



\m -771 1 B{z,t) 



(2) 



where 



A{z,t) = cosh(2772t) + — cosh(277l^)e2*(''2'~''^)^ (3) 
B{z,t) = I^l±^cosh[2(772-r7i)z] (4) 



4771 
{m - mf 



cos [2(772 -^1)^+ cosh [2 (771 +772)2]. 



The two free parameters, 771 and 772, are the imaginary 
parts of the poles in the scattering data, i.e. \j ~ irjj, 
and the residues are related to the poles by the relations 



(j2 ^ Uk=liVj+Vk) 



III. SOLITON SQUEEZING 

After knowing that the general solution for the N = 2 
solitons may have different initial profiles when we change 
the ratio 771/772, we apply the back-propagation method 
[Tsj l to calculate the quantum fluctuations of a full set 
of the N = 2 soliton solutions. To evaluate the quan- 
tum fluctuations around the bound solitons, we replace 
the classical function U{z, t) in Eq. by the quantum- 
fleld operator variable, U{z,t), which satisfles the equal- 
coordinate bosonic commutation relations. Next, we sub- 
stitute the expansion U — Uq + u into Eq. to linearize 
it around the classical solution Uq for the soliton con- 
taining a large number of photons. Then we calculate 
the quantum uncertainty of the output fleld by back- 
propagating the output fleld to the input fleld with the 
assumption that the statistics of the input quantum-fleld 
operators obey the Poisson distribution. In particular, we 
calculate the squeezing ratio, deflned below, of the output 
fleld based on the homodyne detection scheme |26L |27| , 



v^T[{fUt)\u{L,t))] 
^ var[(/L(0|7i(0,t))]^ 



(5) 



where var[-] stands for the variance, fL(t) is the nor- 
malized classical pulse solution in the output with an 
adjustable phase shift. 



hit) = 



Ua{L,t)e 



^C^dt\UoiL,t)\ 



which acts as a local oscillator. The optimal (minimum) 
value of the squeezing ratio R{z) can be chosen by varying 
the parameter 9. When 9 = 0, the in-phase quadrature 
component is detected, and when 9 = 7r/2, the out-of- 
phase quadrature component is detected. 
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FIG. 2; Optimal squeezing ratio vs. propagation distance for 
the N = 2 solitons with different values of the ratio of 7)1/772, 
for 771+772 = 2. Dashed line shows the optimal squeezing ratio 
curve for the case of the A'' = 1 soliton. 



Based on the formulation above, now we calculate the 
optimal quadrature squeezing ratios for a full set of the 
N ^ 2 soliton bound states in Fig. [21 Compared to 
the optimal squeezing ratio of the fundamental soliton 
(dashed line) , all = 2 solitons get more squeezed in the 
beginning of their propagation because they contain more 
energy Moreover, after a certain propagation 

distance, the optimal squeezing ratio of the N — 2 soliton 
state change periodically due to the oscillating behavior 
of the breather, with the period of the N = 2 soliton [23| . 



{v2 + vi){v2 - m) 



(6) 



However, if we compare the optimal squeezing ratios be- 
tween the N — 2 solitons with different ratios rii/r]2, we 
find that solitons with larger ?7i/?72 are more squeezed, 
even all of them have the same energy. 

The reason that a, N = 2 soliton with a larger ?7i/?72 
gets more squeezed can be inducted from the compari- 
son of the optimal squeezing ratio with that of the = 1 
soliton. When the propagation distance is large enough 
(beyond the propagation range shown in the Fig. [2J), 
there is a oscillating tail in the optimal squeezing ratio 
of the iV = 1 soliton coming form the continuum part 
of the noise due to the use of the same pulse profile as 
the local oscillator. But basically the optimal squeezing 
ratio of the fundamental soliton increases monochromati- 
cally along the propagation distance due to the stationary 
characteristic of the pulses. On the contrary, the oscil- 
lation nature of = 2 solitons prevent the increase of 
the optimal squeezing ratio after a certain degree. Since 
a iV = 2 soliton with larger ratio of 771/772 has a longer 
collision period, as shown in Eq. 10, it is this longer 
collision period that makes the pulse to behave more sta- 
tionary, and more squeezed. 
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FIG. 3: Correlation spectra in the frequency domain for A'^ = 
2 soliton bound states, (a)-(c) an initial single- hump soliton, 
771 : 7/2 = 1 : 3, at different propagation distances, 3.1, 3.9, 
and 4.7, respectively, (d)-(f) an initial double-hump soliton, 
rii : ri2 — 1 : 2, ai different propagation distances, 2.3, 3.5, and 
4.7, respectively. Insets show the soliton Fourier components. 
In all cases, 771 -I- 77 = 2. 



In addition, in Fig. |31we present the results of our cal- 
culations of the frequency-domain photon-number corre- 
lation spectra for the N = 2 solitons with two values, 
Vi ■ V2 — ^ '■ 3 and 771 : 772 = 1 : 2. The correlation coef- 
ficients, which are defined through the normally ordered 
covariance. 



AhiAfij :) 



(7) 



are calculated by means of the back-propagation method 
[l5l |. In Eq. O, Ahj is the photon- number fluctuation 
in the i-th slot Acui in the frequency domain, 

Afii — / dt[U {z, u! 

where AJ7(z,a;) is the perturbation of the quantum-field 
operator, U{z,uj) is the classical unperturbed solution, 
and the integral is taken over the given spectral slot. 
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First, we reproduce the results for the photon- number 
correlation spectra of the initial sec/i-like single-hump 
solitons at ?7i : ry2 = 1 : 3, reported earlier by Schmidt et 
al. 0, . A cross-pattern of the anti-correlated com- 
ponents, corresponding to the values C'ij = —1, occurs 
when the solitons merge, as shown in Fig. Ol^b). This 
is the reason why an efhcient number squeezing of the 
NLSE solitons can be produced by spectral filtering that 
removes the noisy spectral components ^2 ■ between, 
the photon-number correlation spectra change periodi- 
cally as the classical soliton profiles Fig. |3Ia-c) It must 
be noted that there are some correlated patterns outside 
the center part of the solitons, even though the amplitude 
of the Fourier components there almost vanishing. These 
correlated components in the far fringes come from the 
breathering dynamics of the N = 2 solitons. 

Now we turn to the case of a initially double-hump 
profile, 771 : 772 = 1 : 2, as presented in Figs. Old-f). 
As can be seen in Fig. |2Ib) and Fig. I^e), the spectra 
for a double-hump pulse contains the same correlation 
patterns as that for a single-hump soliton when all of 
them merge, at z = 3.9, for 771/772 = 1:3, and z ~ 3.7, 
for 771/772 = 1:2. That is what one expects for both 
of them having the same energy and same profile when 
the two solitons merge. And again, the correlation spec- 
trum for a double-hump soliton returns to the same pat- 
tern after a collision period, as shown in Fig. |SJd) and 
Fig.|3Jf). For double-hump solitons, there are significant 
differences in their correlation spectra patterns, see, e.g.. 
Fig. Ola) and Fig. Old). If we only look at the center 
part of the correlation spectra where the soliton Fourier 
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IV. CONCLUSIONS 

We have demonstrated that a two-soliton bound state 
gets more squeezed when it has a double-hump initial 
profile, and this effect is associated with longer soliton 
collision period. We also study the photon-number cor- 
relation spectra of the N = 2 solitons, which reveal 
the anti-correlated patterns which make the soliton to 
get more squeezed. Since such double-hump two-soliton 
states have been generated in experiment, we do expect 
that our theoretical predictions can be readily verified ex- 
perimentally, by generating strongly squeezed states for 
quantum information process and quantum computation. 
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